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Physical conditions in the Protoplanetary Nebula CRL 618 
derived from observations of vibrationally excited HC3N 

F. Wyrowski 1,2 , P. Schilke 2 , S. Thorwirth 3 , K.M. Menten 2 , G. Winnewisser 3 

ABSTRACT 

We used the Effelsberg 100 m and IRAM 30 m telescopes to observe vibra- 
tionally excited cyanoacetylene (HC3N) in several rotational transitions toward 
the proto-planetary nebula CRL 618. Lines from 9 different vibrationally ex- 
cited states with energies ranging up to 1600 K above ground were detected. 
The lines show P Cygni profiles indicating that the HC3N emission originates 
from an expanding and accelerating molecular envelope. The HC3N rotational 
temperature varies with velocity, peaks at 520 K, 3 kms" 1 blue-shifted from the 
systemic velocity and decreases with higher blueshift of the gas. The column 
density of the absorbing HC3N is 3-6 xlO 17 cm -2 . We modeled spectra based on 
spherical models of the expanding envelope which provide an excellent fit to the 
observations, and discuss the implications of the models. Additionally, lines from 
13 C substituted cyanoacetylene were observed. They can be used to constrain 
the 12 C/ 13 C ratio in this source to 10 ± 2. 

Subject headings: stars: individual (CRL 618) - stars: AGB and post-AGB - 
radio lines: stars - stars: mass loss circumstellar matter - line: profiles 



1. Introduction 

At the end of their lives, stars go through a sequence of red giant phases accompanied 
by intense stellar mass loss. This results in atomic and molecular circumstellar envelopes 
enshrouding the central asymptotic giant branch (AGB) stars and leading to bright infrared 
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sources. Eventually the star evolves into a white dwarf hot enough to ionize the enve- 
lope, starting the formation of a planetary nebula (PN). The transition phase between the 
AGB and PN stages is very short; consequently objects evolving through it, so-called proto- 
planetary nebulae (PPNe), are very rare. Their physical and chemical properties define the 
starting conditions for planetary nebula evolution. Therefore, it is of importance to deter- 
mine these characteristics observationally. Reviews on this subject are given in Kwok (1993) 
and Hrivnak (1997). 

CRL 618 is the archetypical PPN. Its circumstellar, expanding envelope emits in a 
multitude of molecular lines (e. g. Bujarrabal et al. 1988) and contains molecules as complex 
as benzene, which was recently discovered with the Infrared Space Observatory by Cernicharo 
et al. (2001); it shows a bipolar morphology /outflow in the optical, CO, and HCN emission 
(Meixner et al. 1998; Neri et al. 1992); the bipolar outflow lobes were imaged at high angular 
resolution with the Hubble Space Telescope (Trammell 2000); a compact, possibly variable 
Hll region is just emerging around the central star (Kwok & Bignell 1984; Mart in- Pint ado 
et al. 1993). 

To study the physical and chemical properties of the circumstellar envelope of CRL 618, 
we started a project to observe the carbon chain molecules HCN, HC 3 N, and HC 5 N in their 
vibrationally ground and excited states. Carbon chain molecules were chosen since they 
might be related to the formation of carbon dust grains. Another objective was that they 
provide many rotational lines of different vibrational excitation in a relatively small spectral 
range, allowing observations of various lines with equal spatial resolution. The vibrationally 
excited levels are populated by IR radiation, and thus provide the interesting information 
about the IR fields and excitation of the gas (Wyrowski et al. 1999) which would not be 
accessible in the infrared itself due to the high dust extinction. 

In this paper, we present an observational study of vibrationally excited HC3N in 
CRL 618 carried out using the Effelsberg 100 m and the IRAM 30 m telescopes. In the 
next section we briefly describe the observing techniques. In Sect. 3, we summarize our 
main results and in Sect. 4 the data will be interpreted in terms of an expanding, molecular 
circumstellar envelope. A discussion of the HC 3 N analysis comparing our new results with 
complementary molecular studies from the literature is then given in Sect. 5. 
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2. 



Observations 



2.1. Effelsberg 100m 



We observed CRL 618 with the Effelsberg 100m telescope of the Max- Planck- Inst it ut 
firr Radioastronomie on 1998 April 23. The front end was the facility 7 mm HEMT receiver 
tuned to the frequency of the HC 3 N v = 0, J = 5 - 4 transition at 45490.313 MHz. Our 
spectrometer was a 8192 channel autocorrelator which we used with 8 overlapping subunits 
of 80 MHz bandwidth to cover a total of 520 MHz. The resulting spectral resolution was 
0.3 MHz after smoothing the data to increase the signal-to-noise ratio. The beam at the 
frequencies of the HC3N lines was 20". Pointing was checked at roughly hourly intervals by 
means of continuum scans on CRL 618 itself. We found the pointing to be accurate to within 
4". Our absolute calibration was based upon continuum scans through W3(OH), 3C84, and 
3C273, assuming flux densities of 3.5, 6, and 25 Jy, respectively, and should be accurate to 



The observations with the IRAM 30m telescope were conducted on 1998 July 17 and 
18. Using the facility SIS receivers simultaneously at 109.3, 136.6, and 209.4 GHz with 
half power beam widths of 22", 17.6", and 11.5", and system temperatures of 200, 350, and 
500 K, respectively. Our spectrometers were an autocorrelator with 1633 channels and 0.31 
MHz resolution and two filter-banks with 1 MHz resolution and 512 channels. All lines were 
observed in the lower sideband with typical sideband rejections of 0.003, 0.025, and 0.04, at 
109.3, 136.6, and 209.4 GHz, respectively. The wobbling secondary mirror was used with a 
beam throw of 200" and a frequency of 0.25 MHz, resulting in spectra with flat baselines. 
The antenna main beam efficiencies, measured on Mars, were 0.63, 0.51 and 0.45, at 109.3, 
136.6, and 209.4 GHz, respectively. We estimate that our absolute calibration uncertainy is 
lower than 25 %. 



25 %. 



2.2. 



IRAM 30m 



3. 



Results 



3.1. 



The 100m observations near 45 GHz 



The integrated, baseline subtracted HC 3 N J = 5 — 4 spectrum of CRL 618, merging 
the different correlator subunits into one spectrum, is shown in the upper panel of Fig. 1. 
The vibrational ground state line of cyanoacetylene shows a P Cygni profile indicating that 
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the line originates from a hot, expanding circumstellar envelope close to the exciting star. 
This situation is illustrated in Fig. 2. In all the vibrationally excited lines the emission part 
is too weak to be detected and only absorption is seen. We have detected pure rotational 
transitions within all fundamental bending modes 1/7, v 6 , z/ 5 , the v^v-j combination mode 
and the Iv-j and 3uj overtones. The line parameters from Gaussian fits to the lines are 
given in Table 1. The emission and absorption features of lines showing P Cygni profiles 
were fitted separately. The brightness of emission lines is given in main beam brightness 
temperature T MB . We report the strength of absorption features in T MB /T C , with T c as the 
continuum temperature, which is proportional to the optical depth of the lines for Tmb/Tc 
much smaller than unity. Another advantage of the T MB /T C scale is its independence on 
the telescope beam size, since the size of the continuum source is considerably smaller then 
the beam (O.lx 0.4 arcsec, Martin-Pintado et al. 1993), as for the millimeter observations 
described in the next section. 



3.2. The IRAM 30m observations 

The lower three panels of Fig. 1 show the spectra of HC 3 N J = 12 — 11, J = 15 — 14, 
and J = 23 — 22. At the frequencies of J = 12 — 11 and J = 15 — 14, even more lines 
lines show P Cygni profiles. Comparison of, for example, all the v$ = If lines demonstrates 
the continuous transition from absorption to emission of the lines at higher frequencies. 
Vibrational satellites corresponding to six different modes of vibration of HC3N can be 
identified: Besides the ground vibrational state, all fundamental bending modes z/7, z/ 6 , z/ 5 
are detected as in the case of J = 5 — 4. Additionally the stretching vibration z/ 4 and the 
combination mode z/ 7 z/ 4 were detected. At 1.3mm (J = 23 — 22), also the 2z/ 6 overtone is 
detected. The latter three vibrational modes correspond to energies from 1300 to 1600 K 
above ground (see vibrational energy levels in Wyrowski et al. 1999). All the line parameters 
are given in Tables 2, 3 and 4. Only heavily blended lines were omitted. Several rotational 
lines of 13 C isotopomers of HC3N in their bending modes u 7 and 2u 7 are observed as well. For 
all the identifications and analyses we used the improved laboratory frequencies of Thorwirth 
et al. (2000, 2001). Moreover, the 109 GHz spectrum shows the J = 41 —40 ground state line 
of HC 5 N along with its vu = le vibrational satellite. A study of this molecule in CRL 618 
is discussed in Thorwirth (2001). 
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3.3. The continuum emission of CRL 618 

With the MPIfR 100m telescope, we observed drift scans centered on CRL 618 to obtain 
its continuum flux density. With the 30m, we estimated millimeter-wavelength continuum 
flux densities of CRL 618 from the total power offsets of the measured spectra, leading to 
consistent results between different scans. This method is reliable since we used the chopping 
secondary mirror so that instrumental and atmospheric drifts are kept small. The measured 
fluxes densities at all observed frequencies are given in Table 5. In addition, we use the 
flux densities obtained during our studies of HC5N and HCN in CRL 618 (Thorwirth 2001, 
Thorwirth et al. 2002). In Fig. 3 we compare the measured fluxes with continuum fluxes from 
the literature. At millimeter wavelengths flux density estimates at similar frequencies agree 
with eachother within a factor of two, but for longer wavelengths the spread is larger. The 
observations below 50 GHz were done using the VLA, the 100m and the Arecibo telescope. 
There is no trend for interferometrically determined flux densities being smaller, so that we 
don't think the interferometer observations miss flux due to the lack of short spacings. There 
is also no clear trend for long term variability: the flux measurements by Turner & Terzian 
(1984) at 2.4 GHz and by Knapp et al. (1995) at 8.4 GHz were performed in 1981 and 1991, 
respectively, and are in accordance with the recent VLA 5 GHz value by Thorwirth et al. 
(2002). On the other hand, variability on shorter timescales might be the case: Kwok & 
Feldman (1981) reported a brightening of CRL 618 by a factor of 2 over a 2-3 yr interval at 
frequencies around 10 GHz. For the following 5 years then, no significant increase could be 
detected (Mart in- Pint ado et al. 1988). Here we only model the continuum emission for the 
frequency range of 45-210 GHz covered in the present study, assuming constant flux densities 
in the short time interval from April to July 1998. The fit to the continuum flux densities 
in Fig. 3 was obtained by calculating the emerging flux from a homogeneous, circular Hn 
region. The radius was fixed to 0.11" (or 190 AU at a distance of 1.7 kpc, Westbrook et 
al. 1975) to agree with the area of the Hn region obtained by Martin-Pintado et al. (1995). 
Then the electron temperature T e and emission measure EM were chosen to fit the data. 
Our estimate for EM agrees with the result by Martin-Pintado et al. (1988). T e from our 
fit is 50% higher. 

4. Analysis 
4.1. Velocities 

In Fig. 1, the frequencies of all HC3N lines are marked assuming a velocity of the gas of- 
20.5 kms -1 . For the emission lines this is a good assumption but it can be seen easily that all 



- 6- 



the absorption features are blue-shifted. This situation is shown in more detail in Fig. 4 for 
all unblended lines with either mostly emission or absorption: for the emission lines, there is 
little variation of velocity with upper energies of the transitions and their average velocity is 
— 20.6±0.2 km s _1 . For the absorption lines the situation is more complicated. Line velocities 
below 1000 K show a rising blueshift with decreasing energies above ground. Above 1000 K, 
the velocities of the absorption is almost constant with a velocity of —27.2 ± 0.2 kms -1 . 
In the framework of an expanding cooling envelope this can be interpreted as lower energy 
lines absorbing mainly from cooler parts of the envelope and higher energy lines from hotter 
parts closer to the exciting star. A constant velocity for lines higher than 1000 K means, 
that no HC 3 N is found closer to the star, hence marking the radius of the envelope where 
the temperature gets small enough that HC 3 N is more effectively produced than destroyed. 

For comparison, NH 3 line velocities taken from Martin-Pintado & Bachiller (1992) are 
shown in Fig. 4 as well. A rising blueshift with decreasing energies can also be seen in NH 3 
but shifted to lower energies above ground. Addtionally, Martin-Pintado et al. (1992) found 
another absorption component at about —56 kms -1 , which was later interpreted (Martin- 
Pintado et al. 1993, 1995) as being due to postshocked gas from the interaction of the high 
velocity outflow with the expanding AGB envelope. In this picture ammonia absorption 
is found symmetrically at ± 15 kms -1 from the terminal velocity of the unshocked AGB 
envelope at -40 kms -1 . This situation is illustrated in Fig. 5 of Mart in- Pint ado et al. (1995). 
In contrast, the HC 3 N absorption profiles are single peaked, suggesting that HC 3 N mainly 
resides in the unshocked gas. It follows then that the ammonia HC component (Martin- 
Pintado et al. 1992) shown in Fig. 4 is a blending of the redshifted lobe of the postshock 
material with the hot unshocked inner part of the AGB envelope, seen in HC 3 N. 



4.2. The high-velocity outflow 

The lines of the ground vibrational state of HC 3 N show broad lines wings. Hence some 
HC 3 N is found in the 200 kms -1 molecular outflow (Cernicharo et al. 1989; Neri et al. 1992). 
We fitted the line wings by masking all blending lines and the fit results are given in Tables 
1-4. The lines with higher excitation do not show any detectable wings, maybe with the 
exception of the J = 5 — 4, vj = If line. Assuming optically thin emission from the HC 3 N 
line wings, the rotational temperature of the high velocity gas is 60 ±7 K, where we included 
the HC 3 N results of Cernicharo et al. (1989) into the fit. The corresponding HC 3 N column 
density averaged over a 22" beam is 46 ± 9 x 10 12 cm -2 . 
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4.3. Temperatures 

A rough estimate of the cyanoacetylene level populations can be obtained by using the 
lines with dominant emission features and plotting the column densities in the levels against 
energy above ground in a Boltzmann diagram assuming optically thin emission. This leads 
to a temperature of 280 ± 20 K which should be considered as the average over all distances 
from the central star. Since the optical depth of the emission lines is not known a priori, 
this temperature is only an upper limit. 

From Fig. 4, it is clear that absorption lines with velocities of -28 to -27 kms" 1 trace 
the hottest component of the gas on the line of sight in front of the continuum emission. 
Using only these vibrationally excited lines, we estimate a temperature of 520 ± 80 K from 
a Boltzmann plot. Since the measured quantity T MB /T C is approximately 1 — exp(r) (for 
Tc » T ex ), the optical depth is known and cannot simulate a higher temperature. Fig. 3 
shows Tq as a function of frequency and it can be seen that for all frequencies considered in 
this study Tq > 1000 K. The reason for the decrease in absoption with increasing frequency 
is that the emission gets stronger and fills up the absorption part of the line profile. 

In a refinement of this method, we divided all the absorption lines into velocity bins of 
2 kms -1 and estimated the temperature in each bin from Boltzmann plots separately. The 
result is shown in Fig. 5 for HC3N (5-4) and (12-11), using a systemic (stellar) velocity t> sys 
of -24.2 kms -1 (see discussion in the next section). A clear increase of the temperature 
with lower velocities, presumably closer to the central star, is evident from the population of 
the vibrationally excited states in both rotational levels, confirming again the interpretation 
that the lines originate from an expanding cooling envelope. 



4.4. A spherical LTE model for the envelope 

To account for the observed P Cygni profiles, a more detailed approach is necessary. We 
therefore modeled an expanding spherical envelope with radial power laws for temperature, 
density, and velocity: 

T = T {r/r ) aT , n = n (r/r ) a " , v = v (r/r ) Q " (1) 

For the Hn region in the center we used the parameters derived from our continuum obser- 
vations in Sect. 3.3. We then calculated for a given molecular line and velocity the resulting 
intensity distribution assuming LTE conditions (see Appendix B), and from that the emerg- 
ing spectra in our observing beams. Additional free parameters are the systemic (stellar) 
velocity v sys of the envelope, the intrinsic line width Av, and the 12 C/ 13 C ratio. For the 



- 8- 



inner radius r we used the radius of the Hn region. 

Our search for the best fit parameters was guided by results of the previous subsections: 
the velocity at the inner radius v + v sys was chosen to lie close to the observed velocity 
of the high energy absorption lines (Fig. 4) and the temperature at the inner radius has 
to be of the order of the temperature found from the Boltzmann plot of these lines. As a 
more quantitative approach, we fit the temperature distribution shown in Fig. 5 with the 
power law given above and adjusted t> sys by minimizing the difference between fit and data, 
leading to t> sys = —24.2 kms -1 . The resulting T is close to the value estimated from the 
high energy absorption lines. The resulting power law index (3 of -0.5 constrains the ratio 
ar/o: v . The ratio a n /a v is best determined by the total widths of the P Cygni profiles, since 
it mainly governs out to which velocities absorption/emission is still detected. Finally, the 
ratio between emission and absorption peaks depends on a v , which determines the actual 
spatial size of the emission, hence a large a v leads to small sizes and hence less emission in 
the modeled spectra. 

At this point, it is necessary to discuss an important limitation on the accuracy of 
the HC 3 N column density determination: the exact temperature dependence of the HC 3 N 
partition function, including both the rotational and the vibrational part, is unknown. In 
a true LTE approximation, all existing energy states have to be included, which leads for 
HC3N due to its relatively low vibrationally excited states to a sharp increase of the partition 
function for temperatures above about 300 K. Collisional rates for the excitation of the 
vibrationally excited states are mostly unknown and the measurements and estimates that 
exist suggest a fast increase of the resulting critical densities with energies above ground 
(Wyrowski et al. 1999). On the other hand, most of the lines might be populated rather by 
infrared pumping than collisions. To reflect these uncertainties, we determined the column 
density using the two possible extreme cases: a) counting all the states in the partition 
function and b) counting only the states up to the highest observed energy above ground. 

For the fit we used only unblended lines with small errors on the line frequencies. The 
resulting best fit parameters are given in Table 6: The temperature exponent and the inner 
temperature are close to the values found by Meixner et al. (1998) for CIT 6 and CRL 618. It 
is clearly hotter and steeper than the temperature structure of an optically thin cloud heated 
by a central source (T =420 K, o:t=-0.4). This is not unexpected for the dense, dusty inner 
part of the CRL 618 envelope: if the dust is still optically thick in the NIR, the envelope 
heats itself to higher temperatures by absorbing its own thermal emission, leading to higher 
temperatures in the inside and a steeper temperature structure (cf. Sect. 4.5). This power 
law can also be compared with results from modeling the envelope of IRC+ 10216, where 
Mamon et al. (1998) found an inner power law index of -0.72 by fitting the model results of 
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Kwan & Linke (1982). Millar et al. (2000) give a temperature index of -0.79 for IRC+10216 
in excellent agreement with our results toward CRL 618. 

For a spherical wind, the equation of continuity requires a density exponent of —(2+a v ). 
We model a steeper density exponent than this for both forms of the partition function. Since 
the HC 3 N density is simply hh 2 times X(HC 3 N), the HC 3 N abundance, this probably means 
that in the flow HC 3 N is destroyed or converted into more complex molecules. 

Fits to data are shown in Fig. 6. Most lines show excellent agreement between data 
and model, which is very surprising given the simple assumptions used. The remaining 
small differences can probably be attributed to violations of the assumptions, notably local 
thermodynamical equilibrium and spherical symmetry (CRL 618 does have a bipolar flow). 

4.5. Exact modeling of the dusty circumstellar envelope 

The results of our LTE analysis of the envelope in the last section can be compared 
with exact modeling of the envelope using calculations of the radiative transport in dusty 
shells (DIRT: http://dustem.astro.umd.edu, Wolfire et al. 1986; DUSTY, Ivezic et al. 1999). 
For the DUSTY modeling we used a 30000 K black body as radiation source, a mixture of 
Draine & Lee (1984) silicates and graphite with a standard MRN grain size distribution, 
and a dust temperature of 800 K at the inner boundary. The density structure was derived 
within DUSTY from full dynamics calculations. The remaining free parameters are the 
optical depth and the thickness of the envelope. Fig. 7 shows the model results compared 
to measured IRAS and JCMT flux densities (Knapp, Sandell & Robson 1993). The optical 
depth r v determines the overall shape of the SED, whereas the thickness s=r(out)/r(in) only 
modifies the submm part of the SED. We subtracted the free-free contribution originating 
from the HII region (cf. Sect. 3.3) from the submm emission. The SED shape of the DUSTY 
models is independent of the luminosity, which just scales the model fluxes up or down. An 
adequate fit to the data is reached using ry = 100, a thickness s= 3000 and a luminosity of 
2 x 10 4 L . r(in) is the dust sublimation radius which is for this model at 0.11 arcsec, hence 
this is the same radius we used in the last section. The outer radius r(out) of 330 arcsec 
is comparable to the 400 arcsec observed by Speck et al. (2000) with ISO. As seen in the 
Figure, a smaller outer radius can reproduce the submillimeter flux densities. 

While these results are interesting in their own right, for our interpretation of the vibra- 
tionally excited HC 3 N lines only the inner part of the envelope is relevant and especially the 
high optical depths of the envelope. We checked the inner temperature structure (thickness 
s<10) of DUSTY and DIRT models with the high optical depth to fit the SED and found 
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that the temperature falls off considerably faster than in the optical thin case. While the 
exact structure of the models differ they can roughly be approximated with a T =-0.8 power 
laws, hence consistent with the results of the last section. 

Since DUSTY allows us also to model the envelope expansion by assuming that it is 
driven by radiation pressure on the dust grains (Ivezic & Elitzur 1995, Elitzur & Ivezic 2001), 
we can compare the velocity structure predicted by DUSTY with our results from the last 
section. The model used to fit the SED leads to a terminal wind velocity of 4 km/s for 
a gas-to-dust ratio r g( j = 200 and a dust grain bulk density p s = 3 g/cm 3 . scales with 
(rgdPs) -1 / 2 (see Elitzur & Ivezic 2001 for details). Hence, explaining a Voo of order 20 km/s, 
as implied by our data, would require very unlikely values for r g d and/or p s . But even then, 
the P Cygni profiles predicted by a DUSTY velocity structure would be very broad and less 
peaked than the observed ones. The reason is that for radiatively driven, optically thick 
winds most of the acceleration takes place in a very thin part of the shell: 90% of is 
reached already at about 2 r in . We therefore conclude that the expansion of the CRL 618 
envelope cannot exclusively be due to radiation pressure and other processes have to be 
invoked, e. g. processes related to the powerful outflow associated with CRL 618. 



4.6. The total absorbing column density and the 12 C/ 13 C ratio 

The observed absorption lines can be used to estimate the HC 3 N column density on 
the line of sight. As discussed in the last section, the major uncertainty here is the par- 
tition function. We use again two extreme cases to constrain the column density: using a 
temperature of 520 K (Sect. 4.3) the HC 3 N column density is 3-6xl0 17 cm~ 2 . 

To estimate the 12 C/ 13 C ratio to a high accuracy, we determined the minimum x 2 -sum 
between data and model for the HC 3 N v 5 = If and the HC 13 CCN J = 15 - 14 v 7 = If 
lines. These lines were chosen, because they appear to be unblended and were observed 
simultaneously with the same receiver. Hence, in the ratio, all calibration uncertainties 
cancel out. Appearing mostly in absorption, these 3mm lines are much more sensitive to 
column density than the corresponding 12 C, 13 C lines at the higher frequencies. The errors in 
the resulting densities were determined by varying the densities around the best fit values. 
This leads to a 12 C/ 13 C ratio of 10 ± 2 and the model fits to all of the HC 13 CCN lines are 
shown in Fig. 8. 
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5. Discussion and Conclusions 

Our results include the best determination of the temperature structure in CRL 618 
so far. The method we used offers several advantages: (i) highly excited absorption lines 
of HC3N in vibrationally excited states are able to probe the hottest, inner part of the 
molecular envelope in a pencil beam. Since vibrationally excited lines with the same J u 
are observed simultaneously, their relative calibration is excellent and the peak temperature 
found is 520 ± 80 K. (ii) Due to the expansion of the envelope, the temperature structure 
can be determined as a function of velocity, (iii) From the emission part of the observed P 
Cygni spectra the spatial extend of the envelope at a certain velocity can be determined, 
which leads then to temperature structure as a function of radius (T oc r -0 - 8 ). 

How does our result compare with previous measurements? Cernicharo et al. (2001) find 
a temperature of 200 K from their ISO mid-infrared observations of HCN and C2H2. They 
were not able to resolve the velocity structure and therefore probe with their absorption 
lines an average over the whole line of sight, hence a lower temperature. In the analysis of 
far-infrared ISO observations of CRL 618 Herpin & Cernicharo (2000) distinguish between 
different shells close to the central star, with temperatures of 1000, 800 and 250 K, decreasing 
away from the star with radii of 0.3, 0.5 and 0.8 arcsec, respectively. The shells with the 
higher temperatures than our estimate seem to probe a photon dominated region, where 
HC 3 N might not be able to exist yet. 

Exact modeling of the dusty circumstellar envelope using DUSTY was able to reproduce 
the temperature structure and the spectral energy disribution, but failed to predict the 
observed velocity structure of the envelope, indicating that radiation pressure alone is not 
sufficient to drive the expansion. 

Isotopic abundance ratios in AGB stars deliver, in principle, valuable information on 
nuclear processes and dredge-up convective processes that bring up enriched material into 
the stars' outer layers from where they finally end up in the circumstellar envelopes (see, 
e.g. Balser, McMullin & Wilson 2002 and references therein). The HCCCN/HC 13 CCN 
abundance ratio of 10 we determine is much lower than the solar system value of 89. Our 
ratio is very similar to the 12 CO/ 13 CO ratio of 11.6 derived by Wannier & Sahai (1987) for 
CRL 618, which is the lowest they found for any of the ten C and one S stars they investigated. 
The 12 CO/ 13 CO ratios for the other stars in their sample range from 14 to 48. In contrast, 
Kahane et al. (1992) derive 12 CO/ 13 CO= 18 and 3.2 for the 1-0 and 2-1 rotational lines, 
respectively, for CRL 618. However, they adopt a lower limit of 30 for [ 12 C/ 13 C], which they 
derive from the [ 12 C 34 S]/[ 13 C 32 S] double ratio. From carbon monoxide observations, Balser 
et al. (2002) determined the 12 C/ 13 C ratio for 11 planetary and protoplanetary nebulae and 
found values between 4 and 32, with a lower limit of 4.6 for CRL 618 being at the lower end. 
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There are thus a number of uncertainties in the determination of the [ 12 C/ 13 C] in CRL 
618. However, observations of C-bearing molecules other than carbon monoxide should be 
useful to get a clearer picture since they might be subjected differently to fractionation and 
photodissociation than the CO isotopomers, where, according to Kahane et al. (1992) these 
processes cancel each other. 

The present study was supported by the Deutsche Forschungsgemeinschaft (DFG) via 
Grand SFB 494. FW was partly supported by the National Science Foundation under Grant 
No. 96-13716. 

A. The HC 3 N Partition Function 

The partition function can be written as the product of its vibrational and rotational 

part: 

Q = Qv-Qr (Al) 

In cold gas the contribution of vibrationally excited states can be neglected and for linear 
molecules the partition function is simply Q = Q r = 1/3 + kT/hB. In hotter environments 
molecules are excited into vibrational states and Q v is given in LTE as 

Q. = £ . . . £ exp ]> + 1) (A2) 

vi=Q V7=0 ^ ' i=5 

where the sums can be evaluated analytically to yield: 

«.=n(i-« (A3) 

(Herzberg 1945). Here i>j is the vibrational quantum number and di denotes the degree of 
degeneracy of vibration % with vibrational frequency of 

For non-LTE conditions, the higher levels might not be populated for two reasons: 1.) 
critical densities for the higher levels are generally larger so that the density might be too 
small to populate the level. 2.) Vibrationally excited states might be excited by infrared 
radiation fields. Higher levels correspond to shorter wavelengths in the NIR, where the radia- 
tion might be too weak to pump the levels. Hence, to estimate a lower limit for the partition 
function we computed Eq. A2 only up to the highest energy observed. The resulting lower 
and upper limits for the vibrational part of the partition function as a function of temper- 
ature are shown in Fig. 9, together with values from laboratory measurements published in 



the Cologne Database for Molecular Spectroscopy (CDMS, http://www.cdms.de/, Miiller et 
al. 2001). 



B. The Radiative Transfer 



To calculate the intensity J tot at a certain position R and velocity v, we integrate over 
the line of sight z at R: 

Itot(R) = J dz [I(R, z) exp(-dr(R, z)) + S(R, z) (1 - exp(-dr(z)))] (Bl) 

with 

S(R, z) = %- 7 . (B2) 

Here, I(R, z) is the intensity integrated up to a position z on the line of sight. S is the source 
function at a position (R, z). T ex is either the kinetic temperature Tk of the molecular gas or 
the electron temperature T e of the embedded HII region. The optical depths at any position 
(R : z) is given as: 



c 2 



dT (R, *) = ^2 An ^ ( X - exp(/iz//T K ) 0o (B3) 



87TZ/ 

with 

, _ 2c^2) ( (v-v(R,z)) 



and the population of the upper states in LTE: 



exp -4 log(2) ^ ( B4 ) 



n u = exp ( ^ ) . (B5) 



Q \kT J 

The optical depth through the HII region is: 

r c = 0.0824 T" 1 - 35 v~ 2A EM (B6) 
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Table 1: HC3N J = 5 — 4 line parameters. The brightness of emission lines is given in Tmb, 
whereas absorption features are given in T MB /T C . Heavily blended lines are omitted. 



Isotopomer 


(v4,v 5 ,v 6 ,v 7 ) 


I 


V 


Tmb or Tmb/Tq 


^LSR 


Av 










(MHz) 


(K) or () 


(kms -1 ) 


(kme 




HC 3 N 


(0,0,0,0) 




45490.314e 


-0.170 (20) 


-38.7(13) 


11.4 


(24) 










0.290 (50) 


-16. 6(6) 


14.1 


(12) 




(0,1,0,0) 


le 


45494. 714e 


-0.120 (20) 


-27.2(17) 


14.0 


(41) 




(0,1,0,0) 


1/ 


45520.454e 


-0.064(14) 


-28.1 (7) 


6.8 


(15) 




(0,0,1,0) 


le 


45564.964e 


-0.130(10) 


-27.3 (3) 


3.5 


(9) 




(0,0,1,0) 


If 


45600.785e 


-0.115(16) 


-27.1 (4) 


2.7 


(5) 




(0,0,0,1) 


le 


45602.171e 


-0.180(16) 


-31.8(4) 


12.5 


(11) 




(0,0,0,1) 


If 


45667.550e 


-0.240 (40) 


-30.8 (4) 


10.2 


(10) 




(0,0,1,1) 





45727.452e 


-0.040 (20) 


-27.7(7) 


6.6 


(11) 




(0,0,1,1) 


2 


45729.306e 












(0,0,0,2) 


0,2 


45779.054e 


-0.080 (30) 


-28.9(2) 


7.3 


(3) 




(0,0,0,3) 


le 


45856.001e 


-0.096(18) 


-27.5 (4) 


3.7 


(7) 



Table 2: HC3N J = 12 — 11 line parameters. (0, 0, 0, 0) wng acounts for the line wings of the 
ground state transitions originating from the high velocity molecular outflow. 

Isotopomer (u 4 , v 5 , v 6 , v 7 ) i Frequency T M b or T M b/7c ^lsr At> 









(MHz) 


(K) or () 


(kms 1 ) 


(kms 


- 1 ) 


HC3N 


(1,0,0,0) 




109023.305 


-0.220 (30) 


-27.6(2) 


3.7 


(3) 


HC 13 CCN 


(0,0,0,1) 


If 


109125.744e 


-0.120 (30) 


-28.5(4) 


3.7 


(9) 


HCC 13 CN 


(0,0,0,1) 


1/ 


109139.707e 


-0.150 (30) 


-29.0(4) 


6.2 


(8) 


HC 3 N 


(0,0,0,0) 




109173.634e 


0.860 (20) 


-20.4(1) 


17.6 


U) 




(0,0,0,0) wn9 




109173.634e 


0.100(10) 


-14.0(14) 


124.0 


(50) 




(0,1,0,0) 


1/ 


109244.222 


-0.330 (50) 


-28.9(2) 


5.6 


(4) 










0.040 (10) 


-19.6(8) 


7.9 


(19) 




(1,0,0,1) 


le 


109306.704e 


-0.150 (30) 


-27.1(4) 


2.7 


(9) 




(0,0,1,0) 


le 


109352.781e 


-0.420 (50) 


-30.4(3) 


6.0 


(5) 










0.090(10) 


-22.2(6) 


9.8 


(14) 




(1,0,0,1) 


1/ 


109469.409e 


-0.100 (30) 


-27.5(4) 


4.1 


(8) 
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Table 3: Same as Table 2 for HC 3 N J = 15 - 14 











1 MB ui J- MB / 1 C 


^LSR 


An 








1 lVJ.±±Zi 1 


l-^J or y 


(kms 


- 1 ) 


( Irvn c 

1 Kill D 1 


HC 13 CCN 


(0,0,0,1) 


1/ 


136404.396e 


-0.166(21) 


-28.9 


(8) 


4.9(12) 










0.060(10) 


-19.8 


(7) 


5.6(12) 


HCC 13 CN 


(0,0,0,1) 


1/ 


136421.848e 


-0.094(25) 


-29.7 


(22) 


5.6(34) 










0.050(10) 


-19.7 


(10) 


7.2(16) 


HC 3 N 


(0,0,0,0) 




136464.411e 


0.900 (20) 


-20.6 


(1) 


16.5(3) 




(0,0,0,0) wns 




136464.411e 


0.200(20) 


-21.9 


(10) 


97.6 (33) 




(0,1,0,0) 


1/ 


136551.798 


-0.306(39) 


-29.5 


(5) 


5.3(9) 










0.100(10) 


-20.3 


(8) 


5.9(16) 




(0,0,1,0) 


le 


136688.252e 


-0.336 (45) 


-31.3 


(10) 


6.2(16) 










0.140(20) 


-19.0 


(8) 


8.8(13) 



Table 4: Same as Table 2 for HC 3 N J = 23 - 22 



Isotopomer 


{v A ,V 5 ,Vq,V 7 ) 


i 


Frequency 


Tmb or Tmb/Tc 


VLSR 


Av 








(MHz) 


(K) or () 


(kms -1 ) 


(kms -1 ) 


HC 13 CCN 


(0,0,0,1) 


If 


209137.396 


0.110(30) 


-20.6(8) 


10.1 (16) 


HCC 13 CN 


(0,0,0,1) 


1/ 


209164.141 


0.130(30) 


-21.6(7) 


9.1 (12) 


HC 3 N 


(0,0,0,0) 




209230.234 


1.120(50) 


-20.8(1) 


17.6(4) 




(0, 0, 0, 0)^^ 




209230.234 


0.370(50) 


-20.6(10) 


94.6 (33) 




(0,1,0,0) 


If 


209362.113 


0.260 (40) 


-20.3(4) 


11.3(9) 




(0,0,1,0) 


le 


209573.178 


0.320 (40) 


-20.9(4) 


10.2(8) 


HC 13 CCN 


(0,0,0,2) 





209597.019° 


0.040 (20) 


-20(2) 


9(2) 




(0,0,0,2) 


2/ 


209628.074 a 










(0,0,0,2) 


2e 


20966 1.546 a 








HCC 13 CN 


(0,0,0,2) 





209634.224° 










(0,0,0,2) 


V 


209663.206 a 










(0,0,0,2) 


2e 


209694.908 a 









a Line parameters were kept the same for all lines 
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Table 5: Observed continuum flux densities of CRL 618. 



Frequency 


Flux 


Error 


(GHz) 


(mJy) 


(mJy) 


4.5 b 


26 


3 


20.2 fe 


425 


60 


24.7 6 


490 


90 


34.8 6 


680 


140 


34.8 C 


740 


140 


40.8 ft 


750 


75 


45. 6 a 


800 


200 


47.2 6 


840 


85 


82.3 C 


1650 


400 


109.2 a 


1800 


200 


136.5" 


1650 


300 


209.0 a 


1350 


340 



a this study 

b Thorwirth et al. (2002) 
c Thorwirth (2001) 



Table 6: Best-fit parameters for expanding envelope model. Note, that the variation in the 
HC3N density parameters reflect the uncertainty in the partition function calculation. 



Description 


Parameter 


Value 


Temperature 


To 


560 K 




ay 


-0.8 


Density 


no 


500 - 1000 cm~ 3 




Ot n 


-4 --5 


Velocity 


Vsys 


-24.2 kins" 1 




VQ 


3.3 kins" 1 




a v 


1.6 




Av 


4 kins -1 


12 C/ 13 C 


X 


10 


Hll region 


T e 


20900 K 




EM 


4.2 x 10 10 cm- 6 pc 




ro 


0.11" 
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4.55 10 



4.56 10 4.57 10 

Rest Frequency (MHz) 



4.58 10 



0.5 - 




J=1 2-1 1 





J_ 



1.09 10 



1.092 
Rest 



10 D 

Frequency (MHz) 



1.094 10 



III 

~_ L „ i 


i | i i i i | i i 
S 3 J=15-14 

o" o" P 

^ ^ °" 
o" o o" 


II 


III 


o" o" 
q d 
o o" 


<-{ ddcboo" 
o~ qqcoqq 


q o 
op o 
oo _ 


z z 
o o 

o o I 
— ^ ^ 




■ zz zzz z 
: uocxjoo 

ou ^ 
: 2 " " uo o 

oocooo 

XX XXX X 




n A .jJ^ ° 

o 

- , i , , r 


I I I I I I 


I I I I I 


, I , , " 



1.364 10 



1.365 10 1.366 10 1.367 10 

Rest Frequency (MHz) 



1.368 10 




j i i 



_J I L 



2.092 10 



2.094 10 2.096 10" 

Rest Frequency (MHz) 



2.098 10 



Fig. I.— HC 3 N spectra toward CRL 618. 
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Fig. 2. — General sketch illustrating the P Cygni line profile formation in an expanding 
circumstellar envelope. The absorption feature is due to the blueshifted part of the envelope 
in front of the strong continuum from the Hll region. 
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10 2 


10 1 
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J 1 1 lJ 


dust (mJy) 










10 1 


10 2 









Frequency (GHz) 

Fig. 3. — Continuum spectrum of CRL 618. The fit of free-free emission to the flux values 
obtained in this study is shown as a dashed lines. Dust emission is shown as a dotted line 
and the sum of both components as a solid line. Literature flux densities were taken from 
Kwok k Bignell (1984), Turner & Terzian (1984), Mart in- Pint ado et al. (1988), Neri et al. 
(1992), Knapp, Sandell, & Robson (1993), Martin-Pintado et al. (1993), Yamamura et al. 
(1994), Hajian, Phillips, & Terzian (1995), Martin-Pintado et al. (1995), Knapp et al. (1995), 
Meixner et al. (1998). Additionally, the free-free emission is given in a temperature scale 
(see Sect. 4.3). 
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Fig. 4. — Velocities from the Gaussian fits to the HC3N lines as a function of the upper 
energies of the transitions. Data from the 1*5 = l/i>7 = 3 Fermi resonance system are at 
about 1000 K. 
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10 J - 



10 2 h 



® HC 3 N (12-11) 
. HC 3 N (5-4) 



T=T x ((v-v ) /v ) a 



T = 561 ± 80 K 
a= -0.51 ± 0.14 



(v-v sys ) (km /s) 



10' 



Fig. 5. — Temperatures estimates from Boltzmann plots of absorption detected in differ- 
ent velocity bins. The temperatures have been determined for HC 3 N (5-4) and (12-11) 
separately to test the reliability of the method. The solid shows a fit for the power law 
temperature distribution discussed in Sect. 4.4. 
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Fig. 6. — Comparison between observations (histogram) and line profiles from the model 
described in Sect. 4.4. The velocity scale is relative to the systemic velocity. 
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■pi — i — i — i — i 1 1 1 pi — i — i — i — i 1 1 r 




A (^m) 

Fig. 7.— DUSTY models of the CRL 618 SED compared with IRAS and SCUBA data 
(Knapp, Sandell & Robson 1993, unfilled squares are corrected for free-free flux shown with 
a separate solid line). The solid line shows a model with ry = 100 and thickness s= 3000. 
The dotted line is the same model with s= 1000 and the dashed line a model with ry = 10. 
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Fig. 8. — Comparison between observations (histogram) and line profiles from the model 
described in Sect. 4.4 for lines of HC 13 CCN. 
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Fig. 9. — Vibrational part Q v of the partition function as discussed in Sect. A compared 
with results from CDMS (Muller et al. 2001). 



